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1. RESISTENCIA DE CONTACTO 


Cuando existe un flujo calorífico entre dos cuerpos en contacto se observa una caída 
de temperatura AT. en las proximidades de la unión debido a que la rugosidad de 
ambas superficies impide el contacto total. En la unión se distinguen puntos de 
contacto y huecos que quedan generalmente llenos de aire, aunque podrían rellenarse 
de cualquier otra sustancia; e incluso en algunas aplicaciones se hace el vacío. 


Situación ideal ye 


(contacto ES 


q contacto 


Situación real 
(contacto parcial) 


(huecos 


Esta irregularidad de la unión da lugar a una resistencia térmica denominada 
resistencia de contacto (normalmente por unidad de área) ¡28 definida como: 


O bien: 


Por simple análisis dimensional, resulta EE mK/W. 


El calor a través de la superficie de contacto se transfiere por diferentes modos: 
2/18 


e Puntos de contacto: por conducción, que no necesariamente tiene que ser 
unidimensional, pues los planos de contacto puntual no siempre son 
paralelos a la superficie de unión de los cuerpos. 


e Huecos: a través del gas por conducción, o por convección cuando la 
diferencia de temperaturas entre las caras del hueco es elevada, pudiendo 
aparecer incluso radiación si dicha temperatura es suficientemente alta. Se 
recuerda que la diferencia de temperaturas entre las caras de una pared 
cualquiera crece con la potencia transmitida, esto es, T, —T, = Rq, donde 


R representa la resistencia de dicha pared. Un efecto similar ocurre con las 
caras de los huecos, esto es, tanto mayor potencia transmitida a través de 
la pared hace aumentar la diferencia de temperatura entre ambas caras del 
hueco. 


Dada la baja conductividad de los gases, la resistencia térmica de los 
huecos suele ser la dominante en la resistencia de contacto. 


Los factores que tienden a reducir la resistencia de contacto son: 
e Una baja rugosidad de las superficies, por favorecer el contacto. 


e La elevación de la presión de la unión por aumentar los puntos de contacto, 
efecto que se ve significativamente mejorado cuando se alcanza el límite de 
fluencia de alguno de los materiales, lo que da lugar a una reducción 
considerable de los huecos. 


e El relleno de los huecos con alguna sustancia de alta conductividad térmica 
(es muy habitual el uso de siliconas conductivas). 


e La potencia transmitida a través de las paredes que se unen, pues según 
se ha visto esto afecta a la diferencia de temperartura de las caras de los 
huecos. Con el aumento de tal diferencia de temperatura se van poniendo 
en marcha nuevos mecanismos de transferencia en los huecos 
(inicialmente conducción, posteriormente convección y finalmente 
convección y radiación). 


La resistencia de contacto puede ser la resistencia dominante cuando se realizan 
uniones de materiales de alta conductividad, como plata, cobre, oro o aluminio, lo cual 
suele ser habitual en la industria aeronáutica y en electrotecnia. 
Algunas consecuencias de la existencia de la resistencia de contacto son: 
e En sistemas de aislamiento reduce las pérdidas o ganancias de calor por lo 
que resulta favorable, aunque en estas aplicaciones nunca se juega con la 
resistencia de contacto para conseguir el nivel de aislamiento requerido. 


e  Ensistemas de intercambio energético resulta perjudicial pues: 


- si la potencia térmica transferida está fijada, la resistencia de contacto 
aumenta la temperatura máxima de la parte más caliente. 


- si el salto térmico entre focos está determinado, la resistencia de 
contacto supone una merma en la potencia transferida. 


Algunos sistemas que permiten reducir la resistencia térmica de contacto son: 
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e Si no se pretende una unión permanente: mediante grasas conductoras 
(silicona conductiva), láminas de metales blandos (plomo, plata, estaño, 
etc). 


e Cuando se requiere una unión permanente: mediante pegado o con 
soldadura. 


Algunos valores típicos de la resistencia de contacto a presiones moderadas y con los 
acabados usuales se dan en la tabla siguiente: 


1 

Materiales R [mi-K/w] 
Cerámica-cerámica 500-3000 
Cerámica-metal 1500-8500 
Grafito-metal 3000-6000 
Acero inoxidable-acero inoxidable 1700-3700 
Aluminio-aluminio 2200-12000 
Acero inoxidable-aluminio 3000-4500 
Cobre-cobre 10000-25000 
Aluminio rugoso-aluminio (condiciones de vacío) - 150 
Hierro-aluminio 4000-40000 
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2. ESPESOR CRÍTICO DE RECUBRIMIENTO 


Sea una tubería (figura 1) de radios interior y exterior r, y rz respectivamente y 
conductividad k, por cuyo interior circula un fluido a temperatura T; con coeficiente de 
convección h;. La tubería se recubre hasta el radio rz con un material de conductividad 
k,. Externamente hay otro fluido a temperatura T. y coeficiente de convección h,. Se 
pretende analizar la evolución de la potencia calorífica q que se transfiere a través de 
las paredes al ir variando el espesor del recubrimiento. Dicha potencia está dada por: 


Recubrimiento 


sj e : 
h,27r,L 27kL 27k,L h,27r,L 


Figura 1 


Puede observarse cómo en la medida que crece el radio exterior rz del recubrimiento 


In r, /r, 


aumenta la resistencia conductiva , al tempo que disminuye la resistencia 


27 k, L 
convectiva 1/h, 27 r,L. Ante esta situación cabe preguntarse si existirá algún espesor 


de recubrimiento que haga mínima la pérdida de calor. Tomando derivadas de q o de 
R¡,, pues el término T;-T., permanece constante, respecto de rz e igualando a cero: 


dR,  d ln, 1 e: 1 ¿l =ñ 
dr, dr,[(27k,L h,2xr,L| 27L|k,r, hrí 


: k ; Dr aa 
de donde se obtiene que rar, =f, ar al que se denomina radio crítico. La 


0 Cc 
drz e 


segunda derivada de la resistencia en r, vale: 


1 1 2 
= + 
Dll ka 


que al ser positiva indica que la resistencia térmica es mínima para rz=r. y, por tanto, la 
potencia térmica transferida máxima, en contra de lo previsto inicialmente. En términos 
del espesor e de recubrimiento, el espesor crítico se define como: 


d?*R 


1e 


2 
dr, RT 


1 
=——>0 
27 L 


r3=T 


e, =F.—F, 


Cc C 


Por un procedimiento análogo, se comprueba que el radio crítico para una pared 


En e k 
esférica es el doble que en la pared cilíndrica, esto es r. = 2. Estos resultados 
e 
tienen implicaciones distintas según se emplee el recubrimiento con función 
aislamiento o disipación: 
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a. Función aislamiento: tiene por objeto reducir el flujo de potencia con el 
exterior. Tal es el caso del transporte por tuberías de fluidos calientes (flujo de 
calor hacia el exterior) o fríos (flujo hacia el interior). 


q q 


q max 


Qo! 


h ñ r* Ta 
Figura 2.a Figura 2.b 


En la figura 2.a se puede apreciar cómo cuando el radio crítico es superior al 
radio exterior de la tubería (r.>r>) la adición de recubrimiento aumenta el flujo 
de calor si rz<r*, o lo disminuye si rz>r*, respecto al calor qy evacuado en 
ausencia de recubrimiento. r* representa el radio de recubrimiento que iguala a 
la potencia calorífica sin recubrimiento. Evidentemente, la función aislamiento 
exigirá como solución válida rz>r*, pues es la que hace disminuir la potencia 
calorífica respecto a la que tiene lugar inicialmente, q/. En esta zona la adición 
de espesores tanto mayores reducirá tanto más la potencia calorífica 
intercambiada. 


Cuando el radio crítico r., es menor que el radio exterior r, de la tubería (figura 
2.b), la adición de cualquier espesor de recubrimiento siempre disminuirá el 
flujo de calor, tanto más cuanto mayor sea su espesor. 


b. Función disipación: permite reducir la temperatura máxima de un cuerpo que 
disipa una potencia calorífica dada, como ocurre en refrigeración de 
componentes eléctricos y electrónicos. 


Considérese que en la configuración de partida se sustituye la tubería por un 
elemento eléctrico macizo de radio r> (figura 3.a), que disipa una potencia q, y 
cuya temperatura máxima es 7; 


Tee Ti R To 
A > 
Conductor 

Figura 3.b 


Figura 3.a 


Haciendo uso de un modelo de resistencias térmicas (figura 3.b), la 
temperatura máxima del elemento se puede expresar como: 


T, =T, +aqR, 


l 
donde R = darias ES Ae 
2ak,L — h,27r,L 


entre los puntos / y e, la cual debe presentar un mínimo en rz=r¿=k/he. 
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representa la resistencia térmica total 


Rmn 


a Te ro 


Figura 4.a 


Te rr 


Figura 4.b 


Al objeto de reducir la temperatura máxima del componente, las únicas 
soluciones válidas son aquellas en las que la resistencia total Ri. con 
recubrimiento es menor que sin él (R¿), pues de lo contrario la temperatura 
máxima del componente sería superior con recubrimiento. Por tanto, cuando 
roa>r2 (figura 4.a) rz debe cumplir r,<rz<r*, requiriéndose menor espesor de 
recubrimiento en el rango r,<rz<r. que en r¿<rz<r*. Es inmediato que la solución 
óptima está justo en rz=r,, para la cual la temperatura máxima del componente 
alcanza su valor más pequeño, por serlo R;.. 


Para r¿<rz no hay solución posible, pues la adición de cualquier espesor de 
recubrimiento haría aumentar la temperatura máxima del elemento. 


Se analiza a continuación la adecuación geométrica de materiales aislantes y no 
aislantes para el recubrimiento de cuerpos de forma cilíndrica de distinto tamaño con 
funciones aislamiento y disipación. Se considerará que un recubrimiento es adecuado 
desde el punto de vista geométrico cuando el espesor requerido sea como máximo del 
orden del tamaño del cuerpo que se recubre. Para ello nos apoyaremos en los radios 
críticos r¿=k./h¿ mostrados en la tabla siguiente, que contempla diferentes materiales 
(aislantes y no aislantes) con convección natural y forzada en aire: 


Aislante No aislante 
Kk,-0,1 W/m-K k,-10 W/m-K 
Convección natural con aire 
h.-10 W/m?.k 90400 Ed 
Convección forzada con aire 
he-100 W/m?.k Eo Hi 


De los resultados de la tabla se deduce que: 


"= Materiales aislantes: apropiados como recubrimiento en cuerpos de 
pequeño y de gran tamaño (l-1 cm y l-1 m, respectivamente), tanto en 
función aislamiento como disipación. 


= Materiales no aislantes: apropiados como recubrimiento sólo en cuerpos de 
gran tamaño (I-1 m), tanto en función aislamiento como disipación. 


El espesor crítico no es exclusivo de las configuraciones cilíndricas o esféricas, sino de 
todas aquellas en las que la adición de recubrimiento suponga un aumento del área de 


convección exterior: 


7/18 


Existe espesor crítico 


Pared Recubrimiento 


No existe espesor crítico 


Cada configuración susceptible de tener espesor crítico, tiene su propia expresión de 
dicho espesor crítico, pudiéndose hablar de radio crítico sólo en cilindros y esferas. 


Resumen: el recubrimiento crítico ofrece una mejora en la refrigeración de 
componentes con generación interna y llama la atención sobre el cuidado que hay que 
llevar cuando se recubren elementos con función aislamiento, pues la potencia térmica 
puede aumentar en lugar de reducirla. 
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3. CONDUCCIÓN TRANSITORIA CON RESISTENCIA INTERNA DESPRECIABLE 


Hay determinados transitorios de calentamiento o enfriamiento de un cuerpo en los 
que puede considerarse que todo él se encuentra a temperatura uniforme. En estos 
casos se puede emplear un método aproximado para calcular la ley temporal de 
temperatura denominado Modelo de la Resistencia Interna Despreciable (MRID). 


Supóngase una placa plana de conductividad k, espesor L y temperaturas extremas T; 
y T2en contacto con un fluido a temperatura T,.<T2. El calor que atraviesa la placa 
puede determinarse como: 


Figura 1. Efecto del número del número de Biot en la distribución de 
temperatura en una pared plana en régimen permanente. 


T,-T, T,-To 
E 1 


q= 


kA hA 


Si se quiere cuantificar el cociente entre la diferencia de temperatura en el cuerpo 
(AT sotido=T1-T2) y en el fluido (AT quido=T 2-1 .-): 


L 
AT sólido y T, -T, _ kA _ Reonductiva-sólido _ hL = Bi (3.1) 
AT uido T, La 1 Reonvectiva Ksótido 

hA 


siendo Bj; el número de Biot, definido como el cociente entre las resistencias térmicas 
de conducción (en el sólido) y de convección. Bi da idea de la caída de temperatura en 
el sólido (T,-T>) respecto a la que tiene lugar en el fluido (T>2-T..). La figura 1 recoge 
diferentes situaciones según el número de Biot. Puede observarse cómo ATsolido €S 
despreciable frente a ATiidó cuando Bi es muy pequeño, mientras que ocurre lo 
contrario (AT ruido insignificante frente a ATsoig)) cuando Bi es elevado. Con valores 
intermedios de Bi, ni ATsotido Ni ATruiao resultan despreciables, verificándose que AT solido 
AT fuido SÍ Bl=1. 
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Durante un transitorio de calentamiento o enfriamiento de un cuerpo su temperatura 
cambia en el tiempo y de un punto a otro en una cantidad que depende de Bi. La figura 
2 muestra la evolución de la temperatura de una pared que se enfría desde una 
temperatura inicial Ty. Durante todo el proceso está en contacto con un fluido a 
temperatura T.,,<Toy. Se observa cómo si Bi es pequeño se puede considerar que 
dentro del cuerpo la temperatura es prácticamente uniforme y, por tanto, 
T(r,t) = T(t), pudiéndose determinar de forma aproximada la ley temporal de la 
temperatura del cuerpo mediante el denominado Modelo de la Resistencia Interna 
Despreciable (RID), que se estudia más adelante. Cuando Bi no es pequeño (tanto si 
está próximo a la unidad como si es mayor que 1), el cálculo T(r,t) del sólido es más 


complejo, y escapa del ámbito de este capítulo. 


Bizl Bi=1 Biz 1 
T =T(t) T=T(x,0 T=T(x,0) 


Figura 2. Efecto del número del número de Biot en la distribución de 
temperatura durante el enfriamiento de una pared plana en un entorno 
convectivo. 


Modelo RID 
Permite determinar la ley de temperatura en el tiempo T(t) de un cuerpo (en 
calentamiento o enfriamiento) siempre que pueda considerase que internamente dicha 


temperatura es uniforme (no hay diferencia significativa de unos puntos a otros). 
Aplicando el Primer Principio a un sólido que satisface esta condición: 


Qentra + Agenera = sale + W + Oña (3.2) 


siendo: 
W Potencia mecánica desarrollada por el cuerpo. En transferencia de calor este 


término está dado por el trabajo de expansión, pudiendo despreciarse dado el 
pequeño valor del coeficiente de dilatación de los sólidos. 
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OT a de : 
=m — = p Vc—- Ritmo al que varía la energía interna del cuerpo de densidad 
t 


Ot 
p, volumen V, calor específico c y temperatura T. 
Estudiemos el caso típico de un cuerpo que experimenta un enfriamiento desde una 
temperatura inicial Ty en un medio fluido a temperatura T., (<T¿) con un coeficiente de 
convección h y sin fuentes ni sumideros internos de calor. La ecuación (2) toma la 
forma: 


0 


sistema 


0=hA(T=T,)+p ve 


ci: T(t=0)=T0 


donde A representa el área de transferencia de calor por convección. 


As: todo el área exterior. 


A2: excluye el área de la base por estar 
térmicamente aislada. 


aislante 
Realizando el cambio de variable: 
8=T-T., > dT=d0 
la ED queda como: 
do __ A 
dt pve 


ci: 8 (t=0)=To-T..=00 


Integrando: 
0 
¡epa ja 
e) pve pvc 
0 
0 hA 
Mm = 2% 
0, pvc 
y por tanto: 
hA 
T-To.  “pve" 
O 


donde h representa el valor medio del coeficiente de convección entre To y la 
temperatura final del cuerpo. La ecuación anterior puede ponerse en la forma: 
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A (3.3) 


siendo t, = Le la denominada constante de tiempo del proceso, cuyo significado es 


el del tiempo necesario para que la diferencia de temperatura  'T-T., caiga hasta el 
36,8% de la diferencia inicial T¿-T... Gráficamente la ley de temperatura es: 


BOTE 
To -T,, 


pCcL, 


La constante de tiempo también puede expresarse como t, = , donde L. es 


la longitud característica, definida como el cociente entre el volumen V del cuerpo y el 
área A de transferemncia de calor. 


La ecuación (3) también es válida para describir el transitorio de un cuerpo cuando 
además de la convección exista intercambio por radiación con el entorno. Para ello 
será preciso que la superficie radiante evolvente tenga la misma temperatura (T.) que 
el fluido (T..) y que dicha superficie evolvente sea muy grande comparada con el 


cuerpo objeto de estudio. En tal caso h debería tomarse como la suma de los 


coeficientes medios de convección h, y radiación h,, esto es: h=h, +h,. 


A pesar de que el uso del valor medio de hen las ecuaciones anteriores da, en 
general, un buen grado de aproximación en la descripción del transitorio, siempre que 
sea posible deberá plantearse e integrarse una ecuación diferencial que no use estos 
u otros valores medios. En el epígrafe siguiente se presenta un ejemplo de esto (caso 
2). 


VALIDEZ DEL MRID 


LD : : hL = 7 
Cuantitativamente el MRID es aplicable siempre que B, = 3 <0/1 (h coeficiente de 


convección y radiación medio, k conductividad del cuerpo, L¿ longitud característica). 
La definición de longitud característica como L. == permite utilizar el MRID con 


cuerpos de forma cualesquiera, ya sean regulares o irregulares. 


Otras situaciones distintas del transitorio con condición de contorno convectiva, vista 
anteriormente, son: 


Caso 2: 
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= Ausencia de convección (h.=0), esto es, cuerpo en el vacío. 

= Existe intercambio radiativo con una superficie evolvente exterior a temperatura 
T., de mucho mayor tamaño que el cuerpo y supuesto ambos cuerpos grises 
difusos. 


¿0 A, (1* 14) pre E 0 
t 
ci: T(t=0)=T0 


Separando variables: 


A 1 
Evaluando ambas integrales | J Za = 3h z arctan a y ordenando: 
xa 4a x+a  2a a 


t=- pre 33!n pde -1n ate =2| arctg pe - arctg lo (3.4) 
4e A. OTE | |T+T, O e 7 


ecuación que no permite evaluar explícitamente T en función de t, To y Te, ni 
puede ser usada cuando T.=0 (radiación al espacio profundo). 


Para T.=0, la resolución de la integral anterior queda como: 


El número de Biot asociado a este proceso sería: 


BS R conDucción Es 
1 


A z 


RADIACIÓN 
donde h, = solT? me rr) 


Caso 3: 
= Existe convección con coeficiente medio h. con un fluido a temperatura T... El 


Cc 
área convectiva es A.. 
= Existe intercambio radiativo con una superficie evolvente exterior a temperatura 
To, igual a la del fluido (T«=T..), siendo las áreas de transferencia de calor por 
convección A, y radiación A, iguales (A.=A;¡=A). La superficie radiante evolvente 
es mucho más grande que el cuerpo en cuestión. 
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= Existe fuente de radiación puntual (foco, radiación solar, etc) con irradiancia G 
[W/m?], que incide sobre una superficie A”, cuya absortancia a esta radiación 


es a. 
= Existe generación interna de calor q, uniformemente distribuida. 


Tos a GA; 
h. — o AGA] + q V= hi +h ATT) + pre 
*— 
¡0 ci: T(t=0)=T, 
Cambio de variable: 
0=T-T., > dt=d6 

por tanto: 

d0 

——+a0-b=0 

dt 

ci: 0(t=0)=T,-T..=00 
donde a= 24 yb MZ y siendo h=h,+h,. La ED puede resolverse 
pve pve 


sumando las soluciones homogénea y particular o bien eliminando la 
inhomogeneidad mediante la transformación: 


o -=9-? d0=d0' 
a 
y, por tanto: 
Ea a0 =0 
dt 


Separando variables e integrando entre 0 y t: 


0 
ln— = —at 
0 


Deshaciendo el cambio de variable: 
b 
T-T,)-= 
(r=T,)-> 


ARO 


Y finalmente: 
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(T-T., A E 
a e pve 
aGA.+qV 
T E 
lina E hA 


hA 
t 


ss TEST Ñ y j 
que se reduce a la ecuación TS PY“ cuando b=0, esto es, si no existe 


0 00 
fuente de radiación puntual (G=0) ni generación interna (q =0). 


El número de Biot asociado a este proceso sería: 


a (h. + h, JNE 


j k 
donde h, = solT? +T2 fr + T,). Y el | tiempo característico: 


p Ve 
+ 


a 


Caso 4: 
= Existe convección con coeficiente medio h,_con un fluido a temperatura T... A; 


es el área de intercambio por convección. 

"= Existe intercambio radiativo con una superficie evolvente exterior a distinta 
temperatura T. que el fluido convectivo T... A, es el área del cuerpo de 
intercambio radiativo. La superficie radiante evolvente es mucho más grande 
que el cuerpo en cuestión. 

= Ausencia de fuente de radiación puntual (G=0). 

= Sin generación interna (q""=0). 

=  |guales áreas de convección y de radiación: A¿=A,=A. 


Llamando h, = solT 2+T7 tr + T,), la ecuación diferencial se puede expresar como: 


h A(T-T,)+h, ATT) + pvc =0 
t 


cuya solución es: 


h.+h, JA 
e r mE a 
(T, T,)+ 5 En (ue T,) 


debiendo considerar Bi la convección y la radiación: 


B, Pee (h. +E JE 


con un tiempo característico: 
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Cuando se utilizan valores medios de los coeficientes de convección h, y de radiación 


h, , las solución de la ecuación diferencial es sólo aproximada el transitorio del cuerpo, 


pero cuenta con la principal ventaja de que es válido para todo tipo de cuerpos, 
independientemente de su forma, lo que lo hace apropiado para cuerpos con forma 
irregular o compleja. Según esto, cuando no se hace uso de valores medios (ecuación 
4), la solución es exacta. Finalmente, hay que indicar que dado que h.¿ no varía 
considerablemente con la temperatura, una solución que sólo dependa de este valor 
medio, es una buena aproximación, que solo presenta una desviación algo mayor 


cuando se utiliza valor medio de h,, esto es, h, . 


Existen soluciones analíticas exactas de transitorios de conducción para cuerpos de 
geometría sencilla (placas, cilindros y esferas por lo general) y con unas determinadas 
condiciones de contorno. Esto está fuera del programa de la asignatura. 


En el caso general de que el cuerpo estuviese sometido a la acción conjunta de de 
una fuente de radiación puntual, convección, intercambio de radiación con una 
superficie evolvente muy grande e incluso con generación interna de calor (q), la 
ecuación del Primer Principio quedaría como: 


di dT 
aGA, +q V=hoAL(T-T.) +00 Art -To)e po 
ci: T(t=0)=T0 


que constituye una ecuación diferencial ordinaria no homogénea, de primer orden y no 
lineal que no puede integrase analíticamente, aunque sí numéricamente. 
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4. CONDUCCIÓN TRANSITORIA: CONTACTO ENTRE SÓLIDOS SEMI-INFINITOS 


Tr 
Distribución de temperaturas entre sólidos semi-infinitos. 


Sean dos sólidos semi-infinitos, inicialmente a temperaturas uniformes Ta y Ts, que se 
ponen en contacto en el instante t=0. Se puede demostrar que el flujo de calor a través 
de la superficie de unión está dado por la expresión: 


"KATA T,)_Kp(1,—T5) 


¡rat mot 


Dr ta aa)”. kpc)y 
T¡-Tg Kk¿L0gB (koc), 


o bien como a = EA (k conductividad, p densidad, c calor específico): 
pe 


de donde: 


Ta fl2o), +5 lO 


 flkoc), + (oc), 


Puede observarse como la temperatura T; de la unión no depende del tiempo, por lo 
que permanece constante mientras los cuerpos están en contacto, estando T; tanto 
más próxima al cuerpo de mayor kpc. Esto justifica la "sensación de más o menos frío" 
cuando tocamos distintas sustancias a menor temperatura que nuestra temperatura 
superficial (Tpie=34%C y Teuerpos=10%C, por ejemplo): 


a. Tierra seca: kpc-2,8-10%  kg*/s*-k*? 1E21:7 6 
b. Corcho: -1,2-10* Ñ T¡=32,4 *C 
c. Madera: -6,9-10* ú T¡=25,7 %C 
d. Acero: -2,8-10% ' T¡=12,1 %C 
e. Cuerpo humano (=agua): -2,5-10% A 


La sensación de frío aumenta en el sentido: b, c, a, d (b parece el menos frío y d el 
que más). En caso que tocásemos estas sustancias con mayor temperatura que la 
corporal, b parecería el menos caliente y d el que más. En general, el producto kpc es 
elevado en los metales y pequeño en los materiales aislantes, por lo que la sensación 
térmica (de frío o calor) es considerablemente mayor en los primeros. 
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